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If dark matter (DM), which is considered to constitute most of the mass of galaxies, is made of 
supersymmetric (SUSY) particles, the centers of galaxies should emit 7-rays produced by their self- 
annihilation. We present accurate estimates of continuum 7-ray fluxes due to neutralino annihilation 
in the central regions of the Milky Way. We use detailed models of our Galaxy, which satisfy available 
observational data, and include some important physical processes, which were previously neglected. 
Our models predict that spatially extended annihilation signal should be detected at high confidence 
levels by incoming experiments assuming that neutralinos make up most of the DM in the Universe 
and that they annihilate according to current SUSY models. 
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There is an increasing hope that the new generation 
of Imaging Atmospheric (7erenkov Telescopes (IACTs) 
would detect in the very near future the 7-ray signal 
coming from the annihilationproducts of the SUSY DM 
in galaxy halos (e.g., 0, E La Hi)- The success of such 
a detection in competition with other indirect or direct 
experiments including accelerators will solve one of the 
most fundamental questions in Astrophysics and Parti- 
cle Physics: the nature of the dark matter. The light- 
est supersymmetric particle (LSP) has been proposed to 
be a suitable candidate for the non-baryonic cold DM 
(d H, see also for a review). The LSP is stable 
in SUSY models where R-parity is conserved 0, 0, 0| 
and its annihilation cross section and mass has the ap- 
propriate relic densities 0, 0| in the range allowed by 
WMAP, i.e. 0.095 < n CDM h 2 < 0.129 We focus 
in the Minimal Supersymmetric extension of the Stan- 
dard Model of particle physics (MSSM) where the LSP 
is the neutralino (x) ■ New upper limits on the neutralino 
mass (m x ) have been estimated due to the constrains on 
the neutralino relic densities Sl x h? provided by WMAP; 
m x < 1500 GeV based on the MSSM including mini- 
mal supergravity (mSUGRA) [l||l|. A lower limit of 
m x ~ 100 GeV has been set by the accelerators [Tlj . 

The number of neutralino annihilations in galaxy ha- 
los and therefore the expected gamma signal arriving at 
the Earth depends not only on the adopted SUSY model 
but also strongly depends on the DM density Pdm(r)- 
This is why the central region r < 200 pc of the Milky 
Way, where the density is the largest, is the favorite site 
to search for this signal. The expected total number of 
continuum 7-ray photons received per unit time and per 
unit area, from a circular aperture on the sky of width 
at (the resolution of the telescope) observing at a given 
direction \?o relative to the center of the Milky Way can 
be written as: 



YSUS = ^ 7 



2m 2 x 



c/(* ) = / J(^)B(n)dn, 



where the factor YSUS (back-spelled SUSY) depends 
only on the physics of annihilating particles and all 
the astrophysical properties (such as the DM distri- 
bution and geometry considerations) appear only in 
the factor Ufi>o). This factor also accounts for the 
beam smearing, where J(^f) — fi os Pd m ( r )dl, dl = 



±rdr/ y r 2 — c?q sin 'J, is the integral of the line-of-sight 

of the square of the DM density along the direction 
and B(£l)dfl is the Gaussian beam of the telescope: 



B(tt)dQ = exp 
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sinO dOdip. 
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The angles 9 and <p are related with the direction of ob- 
servation ^0 and the line-of-sight angle \& by cos \& = 
cos ^0 cos + sin ^0 sin 6 cos ip. We have assumed spher- 
ical symmetry for the DM particles around the Galactic 
Center and that the observer is located in the Calactic 
equatorial plane at a distance (here 8.0 kpc). 

The factor YSUS/47T represents the isotropic proba- 
bility of 7-ray production per unit of DM density. It 
can be determined for any SUSY model given the neu- 
tralino mass m x , the number of continuum 7-ray pho- 
tons 7V 7 emitted per annihilation, with energy above the 
IACT energy threshold (Eth), and the thermally aver- 
age cross section (av) of the DM particles. We can 
then estimate a YSUS parameter range, given a neu- 
tralino mass interval of 100 — 1500 GeV, and a cross 
section (av) interval of 5 x 10 -27 — 3 x 10~ 26 cm 3 s~ 1 
obtained for a sample of MSSM models computed in 
IJ, |21j with relic densities in agreement with the WMAP 
constrains. The number of continuum gamma photons 
produced per annihilation N 7 is obtained by integrat- 
ing the continuum spectrum given by the decay of tt° 
mesons produced in the fragmentation of quarks. It 
can be well approximated by the eq.(18) in i.e. 
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N y = 5/6 x 3 / 2 - 10/3 x + 5 ^+5/(6^) - 10/3, where 
x = E t h/m x . This gives values of the YSUS parameter 
in the range of 10 -34 to 10~ 30 photons GeV~ 2 cm 3 s~ 1 for 
E th horn 1 to 400 GeV. 

A cuspy DM halo pdm(r) oc r~ Q predicted by the sim- 
ulations of the Cold Dark Matter with the cosmological 
constant (ACDMJ is often assumed for the calculations 

of u{* ) (e.g., HasaEEEHilEllil). 

Cosmological N— body simulations indicate that the dis- 
tribution of DM in relaxed halos varies between two 
shapes: the NFW j2^| density profile p(r) = p /x(l+x) 2 , 
x = r/r s with asymptotic slope a = 1 and the steeper 
Moore et al. Q profile p(r) = p Q /x 1 - 5 {l + x) 1 - 5 , a = 1.5. 
The density of DM also depends on two parameters of the 
approximations: the virial mass M v ; r and the concentra- 
tion C = r v ; r /?' s , where r s is the characteristic radius 
of assumed approximation. For Milky Way-size halos 
the average concentration C = 15 and the lcr-variancc 
is Alog(C) = 0.11. For Moore et al profile we define 
concentration as C moore — Cnfw * 1-72. 

Milky Way mass models with adiabatic compression.— 
The predictions for the DM halos are valid only for halos 
without baryons. When normal gas ("baryons") loses its 
energy through radiative processes, it falls to the central 
region of forming galaxy. As the result of this redistri- 
bution of mass, the gravitational potential in the center 
changes substantially. The dark matter must react to 
this deeper potential by moving closer to the center and 
increasing its density. This increase in the DM density is 
often treated using adiabatic invariants. This is justified 
because there is a limit to the time-scale of changes in the 
mass distribution: changes of the potential at a given ra- 
dius cannot happen faster than the dynamical time-scale 
defined by the mass inside the radius. Adiabatic contrac- 
tion of dark matter in a collapsing protogalaxy was used 
already in 1962 H^. In 1980, Zeldovich et al. used 
it to set constraints of properties of elementary particles 
(annihilating massive neutrinos). They were also the first 
to present analytical expression for adiabatic compression 
(for pure radial orbits) and to make numerical tests to 
confirm that the mechanism works. The present form of 
analytical approximation (circular orbits) was introduced 
in [23. If Mi n (r; n ) is the initial distribution of mass (the 
one predicted by cosmological simulations), then the fi- 
nal (after compression and formation of galaxy) mass 
distribution is given by Mg n (r)r = Mi n (r; n )ri n , where 
Affin = A/dm + -A/bar- This approximation was tested 
in numerical simulations |28l l29j | . The approximation as- 
sumes that orbits are circular and, thus M{r) is the mass 
inside the orbit. This is not true for elongated orbits: 
mass M{r) is smaller than the real mass, which a particle 
"feels" when it travels along elongated trajectory. This 
difference in masses requires a relatively small correction: 
mass M should be replaced by the mass inside time- 
averaged radius of trajectories passing through given ra- 
dius r: Affi n ((r))r = M in ((r in ))r in . Wc find the cor- 
rection using Monte Carlo realizations of trajectories in 
the NFW equilibrium halo and finding the time-averaged 



TABLE I: Models and constraints for the Milky Way Galaxy 





Model A 
INr W 


Model H 
Moore et al. 


Constr. 


Virial mass, 1U M@ 


1.0/ 


1.14 




Virial radius, kpc 




o7n 
I l\J 




Halo concentration C 


11 


12 


10.3-21.5 








(1.5a) 


Disk mass, 1O 1O M 


3.7 


4.0 




Disk scale length, kpc 


3.2 


3.5 


2.5-3.5 


Bulge mass, 1O 9 M 


8.0 


8.0 




Black Hole mass, 1O 6 M 


2.6 


2.6 


2.6 


M(< lOOkpc), lO n M 


6.25 


5.8 


7.5 ±2.5 


Stotal, \Z\ < 1-1 kpc 


65 


70 


71 ±6 


at Rq, M pc -2 








Sbaryon at i? ,Af pC~ 2 


47 


53 


48 ±8 


Kirc at 3 kpc, km/s 


203 


205 


200 ±5 


radii (x) » 1.72x°' 82 /(l 


+ 5ie) - 085 


, x = r/r s . 


This ap- 



proximation predicts smaller contraction in the central 
regions, where individual trajectories are very elongated. 
It gives better fits than the standard approximation when 
compared with realistic cosmological simulations [30|. 

In order to make realistic predictions for annihilation 
rates, we construct two detailed models of the Milky Way 
Galaxy by redoing the full analysis of numerous observa- 
tional data collected in [3l]]. The models are compatible 
with the available observational data for the Milky Way 
and their main parameters are given in Table [I] More 
details on the model ingredients and the existing obser- 
vational constrains can be found in |3l| • Fig. ^ presents 
the distribution of mass and density in the models. While 
all observations were included, some of them are more im- 
portant than others. The solar neighborhood is relatively 
well studied and, thus, provides important observational 
constraints. In Table [I] we present two local parameters: 
the total density of matter inside 1.1 kpc E to tai (obtained 
from kinematics of stars) and the surface density of gas 
and stellar components Sbaryon- Circular velocity V4i rc at 
3 kpc distance from the center provides another crucial 
constraint on models as emphasized in |32| . It is difficult 
to estimate errors of this parameter because of uncertain 
contribution of the galactic bar. We use ±5km/s error, 
which is realistic, but it can be even twice larger. Proba- 
bly the most debated constraint is coming from counts of 
microlensing events in the direction of the galactic bulge. 
Our models are expected to have the optical depth of 
microlensing events t = 1.2 — 1.6 x 10~ 6 and, thus, they 
are compatible with the values of r determined recently 
from the observations r = 1 — 1.5 x 10 -6 |3^|, but are 
excluded if r > 2 x 10~ 6 (see [3lJ for a detailed discussion 
on the bulge optical depth in our models). 

Gamma-ray annihilation observability in the Milky 
Way. — The expected neutralino annihilation gamma 
flux, in units of YSUS/10 -32 , can be computed from Eq. 
1 for the compressed DM density profile provided by our 
Milky Way models as a function of the angular distance 
^0 from the Galactic Center. In Fig. [21 we show the pre- 
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FIG. 1: Density and mass profiles for the Milky Way Galaxy. 
Symbols on right panel show observational constraints as 
taken from Klypin et al. f3l| ■ The full and dashed curves 
labeled "Total" are total mass in NFW and Moore et al. mod- 
els with adiabatic compression. DM mass in the NFW model 
is shown by the thick curve. In the central region most of 
the mass is in baryons. Left panel shows the density. The 
top full curve is the density of baryons. The dashed and full 
curves labeled "DM" are for the Moore et al. and NFW mod- 
els with adiabatic compression. The long-dashed curve is the 
uncompressed NFW profile for comparison. 



FIG. 2: Predicted continuum gamma flux as a function of 
distance from the Galactic Center for Models A and B. 
The thick line shows the flux for the compressed NFW DM 
density profile of the Model A and the thin line for the com- 
pressed Moore et al. profile of the Model B. The flux profile 
for the uncompressed NFW profile is also shown for com- 
parison (dotted line). The dashed line give the minimum 
detectable gamma flux F m i n at 5a level for exposure of 250 
hours and E t h=100 GeV for a typical IACT. The inserted pan- 
nel shows the YSUS/10" 32 dependence with the IACT E th . 
For a given E t h, the shadow region scans all the m x , {av) and 
N 7 intervals (see text). 



dieted fluxes. We also show as a comparison the expected 
flux for the uncompressed NFW density profile. The flux 
profiles were determined for a typical IACT of resolution 
at — 0.1° and solid angle Ail — 10 _5 sr. We have multi- 
plied the flux profiles by a factor of 1.7 quoted by Stoehr 
et al. 4] to account for the presence of substructure in- 
side the Milky Way halo [HIH. 

The success of a detection requires that the mini- 
mum detectable gamma flux F m ; n for an exposure of t 
seconds, given an IACT of effective area A e s, angular 
resolution cr t and threshold energy E t h exceeds a sig- 
nificant number M s of standard deviations (M s a) the 
background noise y/Nb, i.e. F m i n A c st / 'ViV& > M s (see, 
e -g-> El 13) The background counts (N b ) due to elec- 
tronic and hadronic (cosmic protons and helium nuclc- 
ids) cosmic ray showers have been estimated using the 
following expressions Q: JV e = 3x 10~ 2 £? th 2 - 3 t A cS Aft, 
N h = 6.1 x 10~ 3 E^ 1 7 tA cS Aft. As an additional back- 
ground, one has to consider also the contamination due 
to isolated muons which depending on the f.o.v. and al- 
titude location of the telescope may be even the dom- 
inant background at some energy range (the "muon 
wall"). Preliminary studies (3(| show that the muon 
background could be as relevant as the hadronic back- 
ground at E t h ^ 100 GeV while it can be effectively 



rejected at lower E t h- The diffuse galactic and extra- 
galactic gamma radiation are negligible compare to this 
background. Gamma point-like sources within the f.o.v 
can be resolved and taken out a posteriori. The E t h of an 
IACT depends on the zenith angle of observation. The 
Galactic Center is visible at different zenith angles by all 
present IACTs (e.g. CANGAROO-III, H.E.S.S., MAGIC, 
VERITAS), but in the best case an E th of about 100 GeV 
can be achieved. Nevertheless, future planed installations 
may reduce the E t h below 10 GeV. The A e g is also sen- 
sible to the zenith angle of observation, here we choose a 
value of 1 x 10 9 cm 2 . This detectability condition will 
allow us to compute the 5c minimum detectable flux 
Fmin in 250 hours of integration with a typical IACT 
of A eS = 1 x 10 9 cm 2 and E th = 100 GeV (dashed line 
in Fig. [21). At a given distance from the Galactic Center 
only the flux values, for a particular model of the Milky 
Way, greater than F m ; n will be detected. The detection 
will be much harder and may result only in a central spot 
in the case of an IACT with higher Eth , as the YSUS pa- 
rameter declines with E t h (see Fig. E} ■ 

The Milky Way models presented here include adia- 
batic compression and likely will result in a detection of 
the annihilation signal no matter what was the initial 
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(uncompressed) DM density profile. For current experi- 
ments this detection will be successful only for the very 
central regions, less than ~ 0.4° in the case of the Model 
A and close to ~ 1° in the case of the Model B. The 
compressed Moore et al. DM profile will provide a more 
extended gamma flux profile. The uncompressed NFW 
DM profile of the Model A will not be detected even in 
the direction of the Galactic Center. On the other hand, 
even the uncompressed Moore et al. profile of the Model 
B will give a positive detection in the very inner regions 
of the Milky Way. 

The effect of the adiabatic compression included in our 
Milky Way mass models, which was previously ignored, 
is a crucial factor. It should be emphasized that for the 
central ~ 3kpc of the Milky Way, where the baryons 
dominate, it does not make sense to use the dark mat- 
ter profiles provided by cosmological N-body simulations: 
the DM must fall into the deep potential well created by 
the collapsed baryons. Thus, the models presented here 
are not extreme: they are the starting point for realistic 
predictions of the annihilation fluxes. One can envision 
few mechanisms to reduce the effect of the compression. 



Transfer of the angular momentum to the dark matter as 
suggested in |3l| is an option. Yet, recent simulations of 
formation of bars indicate that it is difficult to arrange a 
significant transfer of the angular momentum to the dark 
matter. The DM density in the central few parsec can be 
reduced if the central black hole formed by spiraling and 
merging of two black holes [37j ■ It can also be changed 
(probably reduced) by scattering of DM particles by stars 
in the central 2 pc [38|. If this happens, the flux from the 
central 2 pc can be significantly reduced. Yet, it will only 
decrease the amplitude of the central spike. The signal 
from 0.4° still could be detected because it mostly comes 
from distances 30-50 pc, which are much less affected by 
the uncertain physics around the black hole. 
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